Interactions of Ar^+ and metastable Ar^+ with a Si(lOO) surface 
at velocities near the image acceleration limit 

J. Ducree*, J. Mrogenda, E. Reckels, M. Riither, A. Heinen, Ch. Vitt, M. Venier, 

J. Leuker, and H.J. Andra 
Institut fiir Kernphysik, Westfdlische Wilhelms-Universitdt Miinster, Wilhelm-Klemm-Str. 9, 

D-48I49 Miinster, Germany 

R. Diez Muino 

Laboratoire de Physico-Chimie Theorique, C.N.R.S. et Universite de Bordeaux I, 351 Cours de 

Liberation, 33405 Talence Cedex, France 
(February 9, 2008) 



Abstract 

Auger LMM spectra and preliminary model simulations of Ar^+ and 
metastable Ar^"^ ions interacting with a clean monocrystalline n-doped 
Si(lOO) surface are presented. By varying the experimental parameters, sev- 
eral spectroscopic features have been observed providing valuable information 
for the development of an adequate interaction model. On our apparatus the 
ion beam energy can be lowered to almost mere image charge attraction. High 
data acquisition rates could still be maintained yielding an unprecedented sta- 
tistical quality of the Auger spectra. 
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I. INTRODUCTION 



The interactions of highly charged ions (HCI) with surfaces have attracted the strong 
interest of several research groups in the past, getting a strong boost in the last decade due 
to the increasing availability of high performance HCI ion sources and improvement of other 
experimental equipment. 

In recent years, various technological applications of HCI surface collisions have been 
conceived, in particular for the wide field of microscopic and nanoscopic surface modification. 
In order to foster these efforts, a better understanding of the different stages of the scattering 
process has to be attained. Experimentalists hope to take advantage of the quick charge 
exchange processes and the release of the large amount of potential energy stored in the 
HCI. Unfortunately little consent has been accomplished among researchers on the time 
scales and the location of these processes although a comprehensive series of spectra and 
interpretations has been published on this crucial issue so far. 

According to the classical overbarrier model the neutralization of the HCI sets in 
at a critical distance of typically /2c — 15 A in front of the first bulk layer. depends 
on the target work function W and the initial charge q of the HCI. In the region below 
Rci target band electrons are successively captured into resonant ionic Rydberg states with 
n ~ q\J Rao/W . As soon as more than two electrons have been transferred, the highly excited 
hollow atom starts to relax via autoionization processes yielding low-energy electrons. X-ray 
emission is strongly suppressed for light nuclei. Several studies have been carried out 
showing that the overwhelming fraction of the reflected particles is neutral and suggesting 
that the projectile charge q is already compensated on the incoming path. Nevertheless, it is 
commonly accepted by now |^ that the intra-atomic transition rates involved in the cascade 
are by far too slow to perform a complete relaxation of the neutralized HCI in front of the 
surface. 

Autoionization spectra originating from highly charged ions containing initial inner-shell 
vacancies are characterized by a strong and intense low-energy region and a uniquely shaped 
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high-energy branch which can unambiguously be ascribed to intra-atomic transitions involv- 
ing the inner-shell vacancy. Despite of the low transition rates, certain peak structures can 
even be associated with Auger emission from fully populated shells neighboring the initial 
core configuration. 

In order to clarify the evolution from Rydberg populations to fully occupied lower shells 
and motivated by new experimental findings about large fractions of subsurface emis- 
sion within the autoionization spectra, additional interaction mechanisms have been pos- 
tulated and worked into simulations HOHTll. Also a comparison between Auger spectra for 



the same HCI projectile impinging on different target species 0] and a new theoretical 
approach [jl3[ shed new light on the interaction scenario. It seems that the energetic posi- 
tions of target and projectile electronic states play an important role in all direct inner-shell 
filling mechanisms below the surface. After the HCI has penetrated into the bulk, band 
electrons can shield the HCI core charge and directly feed the lower lying hollow atom states 
while generating a plasmon or an electron hole pair |I^ , p!5[] (so called MCV processes). For 
projectiles with high kinetic energies, electrons can be directly transferred from bulk atom 
levels into inner projectile levels yielding a velocity dependent filling rate 0. 

In an attempt to extract information on particular transition types from the spectra, 
experimentalists have analyzed L- Auger spectra of Ar^'^ ions impinging on tungsten |l6Hl9 



copper pOf and gold [||. These early efforts have been obstructed by the large number of 



initial M-shell configurations that had to be considered in the interpretation of the LMM 
spectra with a few distinctive structures only. In recent years, research activities focused 
on K-Auger spectra of hydrogenlike second row ions C^^, N^+, O^"*", F*+, Ne^"^ [3,plH23 



and Ar^^+ instead. Some clearly pronounced peak regions can be identified in most of 
these spectra and assigned to a comparatively small set of initial L-shell configurations. A 
strong systematic dependence of the relative peak intensities of these KLL spectra on the 
experimental conditions has provided valuable information about the contributing ionic shell 
configurations. 

In this paper, we present several series of L- Auger spectra emitted during the interaction 



of Ar^"*" and metastable "heliumlike" Ar®+ ions impinging with beam energies between 8eV 
and 4.6 keV and different experimental geometries on an n-doped Si(lOO) crystal. For the 
first time, we discovered significant modifications of the shape of the autoionization spectra 
for different projectile energies well below 1 keV and for different observation and interaction 
geometries. They include two geometries that largely suppress the detection of all subsurface 
electron emission. The so obtained spectra exhibit a unique peak profile that largely deviates 
from spectra taken under all other experimental geometries. 

These effects are very surprising because in the energy regime below 1 keV all collisional 
M-shell sidefeeding can generally be ruled out and MCV rates can be treated in a static 
approximation. In order to understand the behavior of the spectra at different incident 
energies we developed an interaction model taking into account the special role of the 3d 
subshell, which mediates an efficient M-shell filling via valence band electrons within the 
bulk. Incorporating this model into a Monte Carlo simulation, the observed alterations in 
the subpeak intensities and positions can quahtatively be reproduced. This model is exper- 
imentally supported by a series of L-Auger spectra emitted by metastable Ar*+ projectiles. 
Under the same experimental conditions, the Ar^+ LMM Auger peak structures turn out 
to be amazingly similar to the Ar^+ LMM structures. 

Section || will introduce to the experimental setup implemented for our measurements. 
In Section |T|, we display several sets of autoionization spectra as obtained under specified 
experimental conditions. Section |I^ will describe how the LMM subpeaks in our Ar^+ 
spectra can be assigned to particular groups of intra-atomic transitions. The next Sec- 
tion outlines the basic ingredients of the subsurface interaction model which we employ 



for the simulation of the Auger spectra. In Section 0, we extract information about the 
evolution of the projectile neutralization and Auger emission from the combined analysis of 
experimental observations and the simulation results. Further experimental proof for the 
portrayed interaction mechanism will be given with the discussion of the L-Auger spectra 
of metastable Ar^+ projectiles in |Vli| . Finally, in |V111| , we summarize the basic findings of 
this paper and give a short outlook on future research. 



II. EXPERIMENTAL SETUP 



Highly charged ions are extracted by a fixed voltage of -20kV from an ECR ion source, 
developed in our laboratory. The metallic vacuum chamber of the source can be floated 
on selectable potentials Uq with respect to earth potential. These ions are q/m separated 
by a double focusing sector magnet system including an aberration correction lens. Two 
electrostatic Einzel lenses convey the beam through the intermediate stages of a differentially 
pumped vacuum system which is needed to maintain the pressure gradient between the ECR 
source ~ 1 x lO^^mbar) and the UHV target chamber (p ^ 5 x lO^^^mbar). Before hitting 
the grounded Si wafer, the ions pass through two deceleration lenses which are optimized 
for a maximum of ions deposited on the target surface of approximately Icm^. 

The kinetic ion energy distribution is recorded by an ion spectrometer which is mounted 
on the beam axis close behind the movable target. For Ar^+ and Ar®"*" beams, the full width 
at half maximum never exceeded 2eV per charge. The center of the peak is a measure 
for the kinetic projectile energy after deceleration E-^^^ = qiUq + Up) where Up is the 
plasma potential which builds up between the plasma and the walls of the ECR source. An 
averaged value of Up = 12V has been observed with variations over months of less than ±2V. 
The Si(lOO) surface has been prepared by successive cycles of Ar^+ sputtering at grazing 
incidence and annealing until all impurities have disappeared from AES spectra and good 
LEED patterns have shown up. 

The geometry within the target chamber is displayed in Fig. |l|(a). The beam axis inter- 
sects the target surface at an angle B. Electrons are detected by an electrostatic entrance 
lens followed by a 150° spherical sector analyzer at an angle \1/ with respect to the surface. In 
most measurements we chose + \l/ = 90°. As approaches 0° in Fig. |l](b), the path length 
inside the solid for electrons which are emitted below the surface drastically increases such 
that the detection of above or near surface emission is clearly favored. Due to the chamber 
alignment and the large acceptance angle of rj = 16 ± 6° of our electron spectrometer en- 
trance lens, below-surface emission is always observed, but to a much smaller extent than 
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above or near surface emission. The absolute spectral intensity in the (\E' ~ 0°)-geometries 
greatly diminishes, though. By rotating the target of Fig. |l|(a) around the ion beam axis 
with the surface normal pointing out of the image plane, the condition of + \E' = 90° could 
be relaxed, and geometries with 6 = 5° and \1' = 0° have been achieved. 

The effective incident energy of the ions on the surface is given by -Ej^j^ plus the energy 



gain resulting from the image charge acceleration 



where the work function W equals 4.6eV for our Si target and g = 9. Accordingly, there 
will always remain a minimum incident energy of approx. 29eV leading to an additional 



perpendicular projectile velocity component Aw^ = y2 ■ E-^^/m. Thus the interaction 
period of the ion in front of the surface can principally not be stretched above an upper 
limit depending on q and W even though the original perpendicular velocity component 



Vi_ = \J'^E\ia]_l^ ' cos(9) of the projectile may vanish by selecting Uq = —Up or 6 i-^ 0°. 

When the incident energy E-^^^ is lowered the beam spreads up at the target (Liouville's 
theorem) and incident angles may deviate from their nominal values G. In the energy domain 
-^kin ^ -^im' projectile path is strongly bent by the attractive image acceleration causing 
increased effective incident angles Ogff, especially for small 0. Hence the values given for 
B in this paper are intended to delineate the chamber geometry rather than the effective 
scattering geometry of an individual projectile. 

Projectile penetration depths at the stage of complete neutralization and deexcitation 
can be estimated by multiplying v_[_ by a typical overall interaction time of 10~^^s. With 
E^ = imt>^ expressed in eV, the perpendicular path length Zp^n of the Ar projectile within 



the bulk can be attained from Zpen = 0.22Ax y E^[eV]. This implies that at energies E^ in 
the range of lOOeV, Zp^n stays below one lattice constant amounting to 5.43A for Si. TRIM 
simulations performed for a lOeV and a lOOeV Ar^^ beam impinging on a Si crystal at 
perpendicular incidence yield average lateral ranges of 3±lAand 10±4A, respectively. These 
distances refer to the total penetration depth until the ion is stopped within the bulk. 



III. EXPERIMENTAL OBSERVATIONS 



Using the apparatus described in the preceeding section, we have measured secondary 
electron spectra emitted by Ar^+ and metastable Ar^+ ions during their interaction with 
the Si wafer. In this work we will focus on examining the well defined high-energy L- 
Auger peaks covering the interval between 120eV and 300eV. The spectra also feature a 
low-energy part which extends up to more than lOOeV. The analysis of electron spectra in 
this energy domain is aggravated by a lack of substructures, the superposition of kinetic 
and intra-atomic emission and their sensitivity to stray electromagnetic fields. Regarding 
the high-energy branch, we point out that no background due to kinetic electron emission 
has to be considered for -Ej^jj^ < 121 eV since the collision energies Ecou = -Sj^in + -^im 
are smaller than the lower bound of the spectral region to be examined. By selecting 
Uq = — 20V< Up = 12 ± 2V, we can prevent HCIs from reaching the grounded target. 
Only projectiles that are partially neutralized before the deceleration stages and secondary 
electrons which are generated by collisions of the HCIs with beam transport lens elements 
(these are on negative potentials) can hit the target where they may set free secondary 
electrons. We discovered that both contributions are negligible. 

In Fig. 1^ we present three Ar^"^ spectra measured under B = 45° and with E-^^^ = 9 eV, 
121eV and 1953eV. This and all following spectra are normalized to the total intensity in 
the L-Auger region between 160eV and 240eV. Considering that at maximum one L-Auger 
process per ion takes place, this type of normalization method is suitable to display the 
intensity shifts between L-transition subgroups as discussed in this paper. We note that the 
calibration of the spectra to the absolute beam intensity is prone to errors which emerge 
from the uncertainty in the correction factors compensating geometrical and kinetic effects. 

At first we recognize the general shape of an Ar^+ LMM spectrum featuring a dominant 
peak at 211eV, a broad structure reaching down to about 120eV on the low-energy side and 
a shoulder sitting on the high energy tail of the spectrum. At -Ej^jj^ = 9 eV, this shoulder 
can be resolved into two subpeaks of almost equal height at 224eV and 232eV. Proceeding 
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to higher -E-i^j^, the 232eV-peak disappears and the 224eV-peak gains intensity. Presumably 
due to the poor statistical quality, the latter 232eV-substructure cannot unambiguously be 
identified in de Zwarts |T^ measurement]^ which was taken under the same experimental 



geometry and roughly the same incident energy on a tungsten target. 

The data acquisition statistics of our spectra exhibits a remarkably high quality. Beam 
current shifts during measurements are compensated by an online normalization of the 
spectra to the overall charge current Iq hitting the target. The accumulated counts per leV 
energy channel in the (G = 4:5°,E-^^^ = 121 eV)-spectrum amount to more than 200,000 
at the 211eV-maximum letting the relative error drop below 0.3%. We note that each 
spectrum in Fig. ^ has been recorded in a single five minute run. This is possible due to 
the high current Ig = 125nA on the target which can be converted into a particle current 
Ip by dividing Ig over the projectile charge q and applying a correction factor compensating 
secondary electron emission. Multiplying Ip by an appropriate geometrical factor, it can be 
shown that the overall experimental count rate in the high-energy branch roughly correlates 
to the emission of one high-energy electron per incoming HCl. 

The spectral series of Ar^+ ions impinging on Si (100) with constant Ey^^^ = 121 eV in 
Fig. ^ displays the variation of the relative peak intensities with the experimental geometry. 
Recalling Fig. H, we discover that the presence of a strong 232eV-subpeak is connected to 
minimum perpendicular velocities Vj_. In the measurement under O = 90°, the observation 
angle \E' is very fiat and a second broad peak region evolves around 198eV. Switching to 
the other "grazing observation" alignment at = 5°,\E' = 0°, this structure is preserved 
proving that its presence is related to a small observation angle \l/ rather than the direction 
of incidence 6 or 6gff . 

Under 6 = 5°, the perpendicular projectile penetration into the bulk is principally 



^ There has obviously been a mistake in the calibration of the plot on the energy axis that has 
been corrected in [17|. 



8 



limited to less than one lattice constant. The severe discrepancy between the two spectra 
under = 5° and = 5°, \& = 0° in Fig. ^ illustrates the extreme above-surface sensitivity 
of the (\E' = 0°)-measurements since the physics of interaction is only determined by Ggg 
and -Ej^in which remain constant. We deduce that the broad peak region is generated above 
or at least near the first bulk layer. Because this region looses its weight under \1' = 45° when 
electrons originating from all interaction phases are detected, above-surface processes only 
supply a minor fraction of the total high-energy emission. Nevertheless, the ratio between 
the detected above- and below-surface emission is strongly enhanced at grazing observation 
\1' = 0° and small projectile penetration depths. 

To obtain a quantitative estimate, we ran TRIM calculations [|25| for Ar^+ ions colliding 



with a Si target. The results show that a few percent of the incoming particles are reflected 
for -Ejjjj^ = 121 eV and 2 keV complying with the preceeding interpretation of the (\E' = 0°)- 
spectra. We point out that one has to be careful about adopting these findings for HCI 
beams because the TRIM code solely employs potentials which are strictly speaking only 
valid for onefold ionized ground state projectiles. For incident energies of less than lOeV 
when > -Ej^jj^, the code fails to produce physically meaningful output since it obviously 
misrepresents the potentials evolving from the complex coupling of the HCI-surface system. 
These potentials are decisive for the calculation of the HCI trajectory along the prolonged 
interaction period in front of the surface and the reflection probability. At such low incident 
energies, no experimental data on reflection coefficients of Ar''+ impinging on Si(lOO) are 
available in the literature or refers to grazing incidence conditions where the physics of 
interaction is different despite the similar vertical velocity components. The detection of the 
unique peak profile at grazing observation combined with the oncoming discussion may be 
regarded as indirect experimental evidence for the existence of reflected projectiles. 

The shifts of the upper edge of the 211eV-peak in Fig. ^ can consistently be explained 
by an enhanced below-surface damping of the emitted electrons at O = 90° which is more 
effective than at O = 5°, \1' = 0° due to the higher perpendicular velocity component v±. 

In Fig. ^ we show spectra of Ar^"*" ions impinging on a n-Si(lOO) surface under different 



incident angles with minimal kinetic energies, i.e., = 9 eV. For = 5° and = 45°, 

the spectra are nearly identical reflecting the fact that the self-image attraction is greater 
than the kinetic projectile energy so that the effective angle of incidence Bgg becomes almost 
independent of its original value 6. While approaching perpendicular incidence, the same 
broad region between 160eV and 205eV as in Fig. ^ pops up again. For the two different 
(vp = 0°)-geometries, the main peaks exhibit about the same height. Since v±_ is minimal 
in all four spectra, the upper edge of the 211eV-peak remains sharp and does not shift to 
lower energies due to bulk damping as in Fig. ^. Even more, the high-energy branches above 
211eV coincide almost perfectly. Keeping in mind our particular choice of normalization 
method and the minimum incident energy -Ej^j^ = 9 eV, the latter feature suggests that the 
peak intensity within the high-energy tail region results from above-surface emission which 
is insensitive to bulk damping of the outgoing electrons. 

In Fig. 1^ we present another series of Ar^+ spectra taken at a flxed angle = 90° (i.e., 
\1' = 0°) for different incident energies -Ej^j^- As the point of emissions moves deeper into 
the solid, below-surface contributions are successively flltered out by bulk damping. The 
double-peak proflle transforms into a single unstructured maximum widening to the low- 
energy side as -Ej^j^ increases. The low-energy bounds of the 198eV maximum coincide at 
-^kin = 9 eV and 121eV. The spectrum measured at -Ej^jj^ = 121 eV demonstrates that the 
appearance of the broad peak structure under \1' = 0° and the 232eV-peak occurring solely 
at minimal are obviously not immediately linked to each other. 

The combined analysis of the spectra in Figs. |[]^ renders the following preliminary 
picture which will be supported by further evidence and simulations in the next sections. 
The dominant 211eV-peak originates from below-surface emission since its center moves 
downward, it broadens and its intensity decreases when long path length of the emitted 
electrons through the bulk to the spectrometer entrance can be assumed. Furthermore, 
it does not disappear with growing v±. This also holds for the lower lying part of the 
spectrum. Two equally intense subpeaks on the high-energy shoulder exclusively appear 
when v\_ is minimized. As v±_ increases, the 224eV-peak gains intensity while the 232eV- 
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peak quickly vanishes. This behavior suggests a dependence of the 232eV-intensity on the 
above-surface interaction time even though the resulting emission process may occur after 
surface penetration. 

The broad peak region between 160eV and 205eV under = 0° and Ey^^^^ < 121 eV 
represents near- or above-surface emission since the "detection window" is shallow and the 
chamber geometry favors detection of above-surface transitions at the same time. Subsurface 
contributions are shielded by bulk damping. For reasons that will be given in Section |VII] , 
it is likely that it is made up of a small fraction of above-surface emission from partially 
screened incoming or ionized reflected particles. The preceeding experimental findings will 
play a crucial role in the conception of an interaction model in Section 



IV. ENERGETIC GROUPING OF ATOMIC LMM TRANSITIONS 

In this section we will attribute some spectral features occurring in the energy range 
between 150eV and 300eV to distinct groups of LMM Auger transitions. The energetic 
overlap between neighboring groups will "fortunately" turn out to be sufficiently small such 
that relative peak intensities can be related to the participation of distinguished Auger pro- 
cesses. Furthermore, certain projectile deexcitation mechanisms can definitively be ruled out 
if no intensity is measured in their proper energy range. By merely comparing peak energies, 
we obtain valuable information concerning the HCl-solid interaction which supplement the 
experimental observations of Section |T| before launching any simulation. At the present 
state of research, peak energies can be evaluated more accurately than transition rates for 
the HCI solid system. 



We employ the well known Cowan code |26] in order to simulate configuration energies 



based on spherically symmetrized wave functions for free atoms and ions. In order to 
calculate Auger transition energies within the bulk, we have to take into account the effect of 
the self induced charge cloud consisting of valence band (VB) electrons which surrounds the 
HCI. First approaches have been made on this behalf ]IB|P7| using the density functional 
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theory (DFT). Results show that the nonhnear screening effects due to the electron gas are 
to a good approximation equivalent to the screening by outer shell "spectator" electrons in 
a free atom. 

The hollow atom entering the bulk loses all Rydberg shell electrons due to the screening 
by the target electron gas. The radii of the resonantly populated orbitals are of the order 
of the capture distances, i.e., about lOA and therefore much larger than the Thomas- Fermi 
screening length of less than one Angstr0m as derived in a free electron gas model. Therefore 
all Rydberg levels will be depleted leaving behind the original ls^2s^p^ core configuration 
and possibly some M-, A^- and 0-shell electrons. The target electron gas will swiftly take 
over the role of the outer electrons to screen and so neutralize the HCI charge. A good 
estimate for the reaction rate of the electron gas to the HCI "point charge" perturbation 
is provided by the plasmon frequency which lies in the vicinity of lO^^s^^ for metals. This 
is way above typical rates of the other HCI bulk interaction processes and we can thus 
assume that the HCI core screening by VB electrons is instantaneous. Except for the special 
handling of the transitions with 3d participation, which will be outlined below, all subsurface 
Auger transition energies given in this paper will hence be derived for neutral initial states 
possessing a total amount of q M- and A^-shell electrons and singly ionized final states. 

Let us now look at the grouping of LMM transitions which is plotted in Fig. |^. The 
histogram displays the energetic positions of all LMM transitions originating from initial 
2p^3s^pyd^ configurations {um = x+y+z < 9) of "hollow" Ar^+ atoms which are neutralized 
via q — Hm "spectator" electrons in the A^-shell. Angular momentum coupling as in |^ 
is not taken into account. Each transition is weighted by unity in the plot discarding 
transition rates and statistical factors due to different subshell occupations. For the sake 
of clarity, the whole spectrum is convoluted by a Gaussian function of constant width 2eV. 
This modification evens out conglomerations of Auger lines at certain energies which are an 
artifact of strictly applying the spectator electron approximation. The width is sufficiently 
small not to lead to an additional overlap of LMM subgroup intensities. 

Within the same group. Auger transition energies generally tend to increase steadily 
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with the overall shell population. For comparison, the dotted line in Fig. ^ represents an 
autoionization spectrum of Ar^+ ions impinging on a Si (100) surface at E-^^^ = 121 eV and 
O = \E' = 45° as reproduced from the experimental data in Fig. ^. 

Fig. ^ reveals that LMM Auger transitions involving a free and initially neutral Ar 
atom can cover the energy interval between 166eV (2j9^3s^4s^]9^ i— » 2p^3s°4s^p^) and 267eV 
{2p^3d^ ^ 2p^?)d7). For convenience, the groups of LMM transitions displayed in Fig. ^ 
and the following part of the paper are classified by the angular ^ quantum numbers of two 
participating M-shell electrons. In all cases, the final states are made up of the atomic 
2p level, the remaining M-core states and an appropriate continuum state. For LMM 
processes, we omit the 2p level in our notation. 

The low-energy part of the LMM spectrum can be assigned to 3ss- and 3sp transitions. 
The higher 3sp intensity can be explained by their statistical weight and their 3p contribution 
clearly enhancing the transition rates. The fact that the two small peaks arising in some 
spectra between 190eV and 200eV fall into the 3sp peak region in Fig. ^ might be fortuitous. 
Due to our coarse resolution concerning the energetic grouping, we are not able to ascribe 
these peaks to particular 3sp transitions. 

Several things indicate that the dominant peak region around 211eV is composed of 3pp 
transitions out of a massively occupied M-shell instead of 3s(i transitions the energy range 
of which also covers this peak region. At first, it is intuitively plausible, considering that 
all three bound state wave functions possess the same angular momentum, that the by far 
highest LMM rates are calculated for the "ipp group. Second, the sharp upper edge of the 
2 lleV- maximum resembles the upper boundary of the 3pp curve which is composed of 3pp 
transitions out of a completed M-shell. Due to level filling statistics, a sharp edge is unlikely 
to form if its corresponding transitions take place out of intermediate shell occupations. 
Third, atomic structure calculations yield that 3pp energies accumulate around 211eV for 
all initial 3s^p^# configurations {y + z > 5), regardless of the particular choice of y and 
z. This automatically implies that prior to the majority of all 3pp decays either more than 
seven electrons have to be captured into the M-shell or the induced charge cloud provides 
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an equivalent screening effect. 

According to the LMM grouping in Fig. ^ we can assign the two subpeaks on the 
high-energy shoulder of the LMM-maximum to Ssci- and Spd transitions, respectively. 3pp 
processes are unlikely to contribute to the region above 213eV since they require at least 
one 3s vacancy along with a ninefold occupied M-shell. These initial configurations will 
immediately be converted into 3s^ configuration due to the very fast super Coster-Kronig 
(sCK) decay channel involving three M-shell electron levels. 

The spectral range of the ?)pd peak is cut off at about 235eV and ?)dd transitions do 
obviously not produce enough intensity to appear with a distinct peak region in the spectra. 
These observations provide experimental evidence that the 3(i level cannot be completely 
populated within the bulk and that quick sCK transitions tend to carry 3(i populations into 
lower lying sublevels before LMM transitions take place. The missing structures and the 
spectral range of the high-energy tail extending above 300eV suggest that it consists of the 
large variety of LXY transitions with X,Yg{N, 0} rather than "idd transitions. 

The LMM cut-off at 235eV can be understood by taking a deeper look at the effective 
projectile potential Vgg within the bulk (see Fig. ^ which is deformed with respect to the 
corresponding free ionic Coulomb potential V^l^i- Close to the projectile nucleus r -C aq; the 
effective potential Vgg converges into VqI^i. At intermediate distances r ~ oq, the screening 
of outer levels and the electron gas starts to act on the projectile levels. In this domain V^g 
is well represented by a free atom potential V^'^J^i^'^ which is screened by outer shell spectator 
electrons. All nl subshells with energies E'^' are elevated by a subshell dependent amount of 
Ai?^' with respect to V^^^^. Far away from the nucleus the effective potential V^g merges 
into Vo denoting the bottom of the valence band. 

Fig. 1^ displays the M-sublevel binding energies E]^^ of Ar^"*" as a function of the total 
M-shell population um- The values have been calculated by the Cowan code for spectator 
electron configurations, i.e., for the potential Vq^J^^^. This modeling has proven to yield 
good agreement with experimental and more sophisticated theoretical results in the past. 
In a work by Schippers et al. the main KLL peak energies of the hydrogenlike second 
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row ions C^"*", N^"*", O^^, F^"*" and Ne^"*" have been reproduced. Arnau et al. have 
demonstrated that the spectator electron model complies with DFT calculations including 
nonlinear screening effects for hydrogenlike Ne^+ ions in an Al target. Detailed calculations 
even reveal that the induced charge density tries to mimic the shape of the wavefunctions 
of the neighboring unoccupied atomic level. 

In Fig. I we added 2p binding energies of hydrogenlike C^^ and Ne^"^ as obtained from 
the spectator model and for comparison the DFT calculation for Ne^+ as a function of the 



total L-shell population hl. Following the screening of the atomic spectator electrons 
resembles the screening by the VB electron gas because the inner atomic levels are energet- 
ically separated from the VB much like they are separated from next higher subshell in a 
free atom. This argument holds for the Ar^+ 3 s- and 3p level and also for nearly all L-shell 
levels in hydrogenlike HCIs which are situated between the C^"*" and Ne^^ curves. 

The evolution of the Be? sublevel energies with um in Fig. ^ differs from the lower lying 
subshells, though. We observe that the 3(i level binding energies are significantly closer to 
the VB and grow above Vq as soon as more than five electrons populate the M-shell. We 
performed a DFT calculation showing that "id electrons are already lost to the VB continuum 
for Um > 4. The spectral cut-off in the Spc? transition domain in Fig. ^ can now be explained 
by omitting all contributions from "ipd transitions with tzm > 4. Aiming to correct for the 
shape of V^g which largely deviates from V^o^*;^" for E]^^ ~ Vq (see Figs. |^ and ||), we shift 
the atomic Sc? level to Vq for hm < 4 to attain higher transition energies compared to the 
mere spectator electron model. In this manner we derive the experimental "isd- and "ipd 
peak positions on the high-energy shoulder within an accuracy of 2% and 1%, respectively. 

V. MONTE CARLO SIMULATION OF THE SUBSURFACE INTERACTION 

PHASE 

In order to elucidate the interaction mechanism which eventually generates the measured 
spectra we worked out a Monte Carlo simulation |^ . Our goal was to reproduce the intensity 
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shifts of the observed spectra for different incident energies in Fig. & On the analogy of 



previous simulations by Schippers et al. 0, Page et al. |[T^ and Stolterfoht et al. [|Tl|] on the 
L-shell filling of hydrogenlike highly charged ions at metal surfaces, we only keep track of the 
populations of the two innermost projectile shells containing at least one vacancy and focus 
on the most dominant transition rates. The ionic cores are neutralized by A^-shell spectator 
electrons. Among all intra-atomic Auger processes, only those yielding an electron above 
the vacuum level are considered. 

During the simulation, the three M-subshell populations are recorded continuously. 
Transition rates, transition energies and sublevel energies are evaluated dynamically at each 
iteration step according to the particular {n^sln^pln^a} configuration. From one step to the 
next, only the fastest transition which is derived statistically from its nominal rate takes 
place. The Monte Carlo method implies the averaging of the simulation results over a suffi- 
cient amount of projectiles. We find that the simulated spectra converge after ~ 1 x 10^ 
particle runs and chose A^ = 1 x 10^. In our implementation of the subsurface cascade, each 
particle is started at the first bulk layer with a fixed angle of incidence G = 45° and energy 
-^kin' -^^^ -^kin ~ ^'^^ ^ assume an initially empty M-shell. 

A. Intra- Atomic rates 

The LMM rates are evaluated by a fit expression proposed by Larkins for free 
multiply ionized atoms possessing no A^-shell spectator electrons. Accordingly, if one or two 
of the n electrons of a subshell which could contain hq electrons are involved in an Auger 
process, the Auger rate calculated using the formulae appropriate for a filled shell F^^/^^*^ is 

(2 1 1 1 

reduced by n/no or [n{n — l)]/[no(no — 1)], respectively. Values for F^^^ o^ily supplied 
for 3ss-, 3sp- and 3pp transitions in the literature which account for the greatest part of the 
overall LMM intensity in the literature. For 3sd-, 3pd- and 3dd transitions, we scale the 
LMM rates Fl^jl^"-^ to reproduce the experimental peak heights. Table | lists the six F^^J^^*^ 
rates which are held constant for different simulations. 
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These LMM rates should not be greatly affected by the embedding of the HCI into 
the electron gas because they chiefly depend on the radii of the participating M-subshells 
which remain fairly unchanged. To show this we recall that the shape of the induced charge 
cloud is similar to the A^-shell. Within the hydrogen atom approximation, the radii of the 
screening cloud r^c and the atomic shells (schematically inserted in Fig. both scale with 
(n — 1)^(1 + 1[1 — ^^^^^"^2 ]}- The ratio r^c/rsp = 2.5 with sc = Ap has to be related to the ratio 
^sp/r^s amounting to 0.83. Due to its great extension, the screening electron cloud should 
therefore have a minor impact on the M-shell orbitals and hence on the LMM rates given 
in Table |. 

Since we do not resolve A^-sublevels, Coster-Kronig MMN transitions have to be handled 
by a global base rate for each M-level pair. In a simple approach, we weight each MMN 
base rate by the initial M-sublevel occupation and final state vacancies such that the average 
rate amounts to 3 x lO^'^s"^. For the purposes of this paper, only the order of magnitude with 
respect to the other transition types matters. We remark that Armen and Larkins have 
calculated transition rates for MMN decay channels which are of the order of 4 x lO^^s"^, 
depending strongly on the angular coupling. This is in sufficiently good agreement with 
our assumption. Only MMN transitions with a final state above the continuum level are 
included leaving over solely (3s)(3(i)N transitions which are of particular importance for the 
initial phase of the interaction. 

The sCK MMM rates are known to be 10 to 100 times faster than any rates for Auger 
transitions possessing the initial and final holes in different principal shells. In our simulation, 
they mainly serve to regroup any M-shell configuration into the appropriate M-shell ground 
state before LMM transitions take place. To achieve this, we utilize a base rate of 1 x lO^^s"^ 
which is scaled by the M-subshell occupation statistics. In Table y we put together the 
average number of MMM processes per particle and the average M-sublevel occupation 
at the time of MMM emission for the two sCK transitions which are relevant for our 
simulation. 
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B. MCV filling within the bulk 



Target levels below Vq can be filled by transitions involving electrons of valence band 
states (C) which are perturbed by the ionic core. The energy gain is conveyed either onto 
another VB electron which is emitted into the continuum or a collective excitation (plasmon) 
is created in the medium. The theoretical approach including the charge displacement in 
the description of the excited outgoing electrons is much more complicated and, at present, 
only unperturbed valence band states (V) are included in the calculations |n,T^. The VB 
electrons take on the role of outer shells in a free atom. 

Using DFT to describe the interaction between the ion and the metal valence band and 
following the same scheme as in |]l5l, we have derived MCV rates for the Ar^+-Si system. 



Table |T| lists the rates per spin state T^*-^^ into the three M-sublevels with the number of 
initial M-shell electrons um as parameter. These MCV rates still have to be multiplied by 
the number of unoccupied final states in the particular M-sublevel to attain actual transition 
rates between two atomic configurations. Tf^^"^ denotes the overall MCV rate into the M- 
shell after carrying out the appropriate statistics. Since sCK transitions are much faster than 
MCVs (cf. Table |ID, we only consider "Coster-Kronig final states" as initial configurations 
in the DFT calculation. The transition rates are independent of the projectile velocity Vp 
equaling their static values for all incident energies occurring in this work. 

Table |T| reveals that Tg^*^^ assumes by far the highest values. Taking into account the 
high degeneracy of the Sc? level, effective rates T^^^^ exceed T^^^^ and Tg^*^^ by more than 
one and two orders of magnitude, respectively. With increasing um-, MCV transfer into the 
3p state accelerates reaching the Tg^*^^ values at low um- This is important considering that 
for hm > 4: the 3d shell vanishes and MCVs into the 3p level constitute the most effective 
M-shell filling mechanism which is eventually responsible for the formation of the dominant 
211eV-peak. 
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C. Collisional filling 



For projectile energies above 1 keV, sidefeeding into the HCI M-shell due to direct 
electron transfer from target atom core levels supplies a velocity dependent filling rate. 
The transfer crossection increases with the energetic vicinity of inner projectile and target 
states 1^ which is maximum for the Ar^+ 3s level with the 2p bulk level of Si possessing 
= 109 eV (cf. Fig g). Experimentally, a Si target LMM Auger peak for spectra with 
-^kin — ■'- ""^^ observed which is directly connected to the vacancy transfer. For 2 keV 
projectiles traveling through a silicon crystal in (lOO)-direction, collisional filling supplies a 
3s sidefeeding rate of Fg"" = Vp/d = 1.8 x lO^^s"^ going on the assumption of one electron 
transfer per collision. Within the energy range below 1 keV, collision frequencies are small 
and the distance of closest approach is too large even for head-on collisions to allow a 
sufficient level crossing for sidefeeding pT| . 



D. Simulation of the 121eV- and 2 keV-spectra 

In Fig. 1^ we plot the experimental spectra from Fig. |^ into three subplots and compare 
them with our simulation results which are convoluted by a Gaussian function of 3eV width. 
In this section we look at the Ar^+ spectra and postpone the discussion of the Ar^+ spectra 
which are displayed in the same plot to Section |V11| . The difference between the simulated 
spectra in (a) and (b) stems from collisional filling which is exclusively enabled for E-^^^ = 
2 keV. In addition, we performed a convolution of the 2 keV-spectrum with an exponential 
function with a decay length of 3 a.u. to compensate for elastic and inelastic energy losses 
of electrons on their way through the bulk region. For -Ej^j^ < 2 keV, this damping becomes 
negligible due to the shallow projectile penetration. 

The intensity ratios among the different LMM subgroups and their peak positions are 
approximately reproduced. The 3pp region displays too much intensity, though which might 
be caused by the LMM rate fit formula (cf. Section [V A| ) overestimating the 3pp rates for 
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high M-populations, see also (Table VI). The 3ss intensity is clearly too low suggesting 
that other transitions types not considered in our model may contribute to this region. The 
enhancement of the 3sd peak parallel to the disappearance of the 3pd peak and intensity gain 
of the 3sp region towards the -Ej^jj^ = 2 keV-spectrum as a consequence of the coUisional 
filling can nicely be observed (cf. Table |ip. The average M-sublevel populations at the 
time of LMM emission (cf. Table |) indicate that the high-energy shoulder is generated 
along the early subsurface interaction phase. On the other hand, the dominant 3pp peak 
occurs at high M-populations benefitting from the growing MCV rates into the 3p level 
and the disappearance of the 3d level towards high um- The missing 3dd intensity confirms 
the presence of the fast MMN and MMM decay channels which inhibit the buildup of 3d 
populations larger than one. 

In the experimental spectra, the low-energy tail displays much less structure than the 
simulation indicating that the mere spectator electron model might be incomplete. We 
carried out other simulations where 20% of the LMM transitions start out from singly 
ionized initial configurations such that the peak regions loose part of their intensity to 
the low-energy side. Doing so the intensity dip around 200eV gets partially ironed out 
and the low-energy tail stretches beyond 160eV. A similar effect could be induced by the 
consideration of L2,3MMM double Auger processes |^2| for which Carlson and Krause ||3^ 
measured a relative contribution to all radiationless transitions of 10±2% and energy shifts 
of more than lOeV |0. For the sake of the clarity of the displayed simulation results we 
did not implement this correction in Fig. |^. 



E. Simulation for a statistical initial M-population 

It is very surprising that by reducing the incident energy from about 121eV to 9eV, a 
significant shift in the relative peak intensities still takes place. On the one hand velocity 
dependent below-surface filling can be ruled out in this energy domain, on the other hand this 
effect must originate from different subshell populations at the time of LMM emission. Let 
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us assume for the moment that individual M-subshells of each particle are filled statistically 
(by a Poisson distribution which is cut off at the subshell degeneracy) at the first bulk layer 
according to their respective degeneracy, i.e., (r;,3£)=2/18, 6/18 and 10/18, multiplied by the 
mean total M-shell population {um) for £ = 3s-, 3p- and ?>d level, respectively. In Fig. |^ we 
present results of a Monte Carlo simulation with {um) = 2. 

For a greater part of these initial configurations, new M-shell redistribution channels 
open up via MMNs and sCKs which are energetically forbidden for hm = and carry part 
of the 3d population immediately into the 3p- rather than the 3s level. The simulations 
in Fig. ^(b,c) and Table | indeed reproduce the intensity shift from the 3sd peak to the 
3pd peak at 232eV going from -Ej^j^ =121eV to 9eV. We remark that this simple model 
of an initial M-shell population before bulk penetration does not hold exactly for the Ar*+ 
simulation where we set n^s = 1, (^sp) = 1 and {risd) = 1- We are going to provide a physical 



motivation for the model in Section VII 



VI. THE EVOLUTION OF THE SUBSURFACE CASCADE 



According to the experimental clues and arguments of Sections |T| and the over- 
whelming part of the high-energy branch originates from below-surface emission. For this 
phase, we designed the simulation presented in the previous section. In the following we de- 
scribe the evolution of the subsurface cascade on the basis of the simulation results combined 
with the experimental data. 

As the HCI penetrates into the crystal bulk region all electrons that have previously 
been captured into outer Rydberg levels will be lost and band electrons will neutralize the 
core charge over a distance of roughly the Debye screening length of the electron gas. Thus 
a second generation of hollow atoms emerges within the bulk. 

Prior to any electron capture, the 0-shell of the Ar^+ core is the uppermost ionic shell 
to still fit below Vq. As long as not more than two electrons populate inner levels, solely 
XCV transitions (with Xg{L, M, N, O}) can proceed. Since the XCV transition probability 
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increases with the effective screening and degeneracy of the final level, XCVs preferably 
populate the 0-shell. Before any significant NOO and MNO Auger emission can take place, 
the rapid XCV filling successively pushes the O- and A^-shell above Vq. This period is 
accompanied by LCV, LNO, LMN transitions etc. creating the smoothly decreasing part 
of the spectrum above the 3j9(i edge. We note that this early phase of the neutralization 
may already start before complete bulk penetration when the projectile travels through the 
vacuum tail of the valence band. 

The loss of whole atomic shells into the valence band stops when the M-shell is reached. 
At this point of the scenario, a low M-shell population with a statistical preference for the 
3(i level (due to its high degeneracy) is likely to occur. MMN-CY. processes transfer these 
3(i electrons quickly into the 3s level before a large 3(i population can accumulate. Other 
MMN transitions (3p)(3(i)N and later (3s)(3p)N are energetically forbidden. This M-shell 
redistribution is accelerated by high speed sCK processes with rates of the order of lO^^s"^. 
Whereas ?)sdd transitions are immediately possible, MMM transitions into the 3p level 
require hm > 3. Along this early M-shell redistribution phase the M-population remains 
fairly constant at hm — 2 + n^s, though because one M-electron is lost along each MMM 
process. Thus 3sd-LMM processes out of initial 3s^d constellations are characteristic for 
this phase causing the 224eV-peak in the experimental spectra. It lasts comparatively long 



because the condition um — 2 keeps the MCV rates (cf. Table p| ) minimal. 

The Ar^+ core will always be surrounded by an induced VB charge cloud (C) because 
the number of bound states nf, below Vq is smaller than the projectile core charge q = 9 
(Fig. H). Hence MCV processes continue to populate empty M-levels faster and faster with 
increasing um- As soon as um > 3 is satisfied, 3pdd sCKs become energetically possible 
and a 3p population builds up while the 3d population remains approximately at one due to 
the presence of the MMN and sCK decay channels. 3dd transitions require the transient 
formation of very unstable M-shell configuration that are unlikely to occur so they do 
not appear in the spectra. At um > 4, the 3d level vanishes into the valence band thus 
interrupting further 3sd- and 3pd emission. Since the 3pp-LMM transitions possess much 
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higher rates than any other LMM transitions they clearly prevail during this later stage of 
the subsurface interaction. 

The dominant peak which is centered at 211eV for -Ej^jj^ < 121 eV in Fig. |^ corresponding 
to 3pp transitions with > 7 provides evidence for the described mechanism, in particular 
for the high MCV rates into the 3d- and later the 3p level. The intensity gain of the 3sd 
peak with respect to the 3pd peak for high E-^^^ is consistent with the greater time window 
of the former transition during the early interaction phase. This effect furthermore verifies 
the assumption of collisional sidefeeding into the 3s level and therefore the 224eV-peak 
assignment by itself. All phases are accompanied by 3ss- and 3sp-LM M transitions which 
constitute the low-energy tail and the region around the two faint subpeaks between about 
180eV and 200eV, respectively. 



VII. SPECTRA OF METASTABLE Ar^+ PROJECTILES 

Seeking to extract additional experimental evidence for the described Ar^+ interaction 
mechanism, we performed a series of measurements involving metastable (2p^3s) Ar^+ ions 
colliding at B = 45° and various kinetic energies with a Si crystal (Fig. |10D. A straight 
comparison with the corresponding Ar^^ series in Fig. |] shows that the general shape of the 
spectra is unaffected by the additional 3s electron except for a slight enhancement of the 
3ss- and 3sp intensities. In fact, the only new structure observed is a small peak arising at 
247eV for -Ej^j^ = 8 eV and generally for lowest perpendicular projectile velocities which 
can also be deduced from Fig. 111]. 



The 247eV-peak has been discussed in detail in |^ along with corresponding peaks which 
occur under similar conditions in the spectra of second row ions in ls2s configurations. It can 
be assigned to so called LMVw transitions in the course of which the 3s electron jumps into 
the 2p vacancy. The emitted electron comes from a level possessing a binding energy which 
equals the target work function W = 4.6 eV for silicon. Due to the shape of the subsurface 
potential V^g (Fig. |^, these levels cannot exist after projectile penetration into the bulk 
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occurred. As mentioned earlier in Section |I^, the strong decrease in spectral intensity above 
235eV gives evidence for this assertion. 

The identification of an above-surface LMVw peak suggests that inner atomic shells 
Xg{M, N, O, . . .} could be partially filled before bulk penetration by an autoionization 
process XVvi/Vi^. We mentioned earlier that also MCV set in with continuously increas- 
ing rates as the HCI travels through the vacuum tail of the valence band. Compared to the 
MCV filling within the bulk, these near-surface M-shell filling channels are likely to proceed 
significantly slower, though. Since sCK processes require certain minimum M-shell popu- 
lations they are widely inhibited for these constellations. One can thus expect that a very 
slow projectile might enter the bulk region with a low M-shell population (um) — 2 favoring 
the 3d level due to its degeneracy. This way we can motivate the ansatz for the simulation of 
the spectra at minimum -Ej^j^ in Section even though an explicit experimental evidence 
is still missing. 

The astonishing similarity of the rest of the Ar^+ and the Ar^+ data bears out our 
previous assumption of fast MCV, MMN and sCK processes within the bulk redistributing 
any M-shell population swiftly into the 3s level. In order to compensate for the additional 
3s electron in Ar®+ M-shell, sCKs have to proceed before an LMM transition takes place. 
This automatically implies that the M-shell must be sufficiently populated and quickly 
replenished at this point. Because a large M-shell population far in front of the surface 
would be in contradiction to all previous experiments we can exempt the above-surface zone 
as the origin of the emitted electrons. This obviously also holds for the 3pd peak at 232eV. 

We made use of the great correspondence of the Ar^"*" and the Ar^+ spectra to check 
the mechanisms and rates entering our interaction model. For the simulations on Ar*+ 
projectiles which are also shown in Fig. ^ we kept the same transitions types and rates 
but added an 3s electron to the initial M-shell population. Within the accuracy of our 
interaction model, the similarity of the two series is well reproduced. 
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VIII. SUMMARY AND DISCUSSION 



In this work we have presented detailed experimental results on the interaction of Ar^+ 
and metastable Ar^+ ions impinging on a Si(lOO) crystal. Doing so we focused on autoion- 
ization spectra measured at low impact energies. In this energy domain, we identified several 
new spectral features which alter with the perpendicular projectile velocity component and 
with the angle of incidence and observation. A consistent interaction model has been sug- 
gested for which MCV processes and the energetic vicinity of the Ar^+ 3d subshell to the 
bottom of the silicon valence band play a decisive role. 

The subsurface interaction phase has been simulated using a Monte Carlo code. Feeding 
the code with realistic transition rates, we have been able to reconstruct the experimental 
peak positions and intensity shifts for different projectile energies. Our results give indirect 
evidence for a very effective below-surface MCV filling as postulated by theory. In contrast 
to KLL spectra of hydrogenlike second row ions impinging on metal surfaces, the main 
intensity of the Ar^+ LMM spectra is located on the high-energy side of the peak region 
corresponding to a massively occupied 3p subshell. We demonstrated that this peculiar shape 
of the high-energy region is linked to the special role of the 3d subshell which mediates a 
fast M-shell filling in the beginning and later disappears due to the screening of the valence 
band electron gas. 

We presented spectra measured at small observation angles with respect to the surface 
parallel. They contain a high intensity peak region which most likely originates from Auger 
emission of incoming or reflected projectiles which do not experience the full bulk screening, 
yet. In addition, we spotted a distinct peak in the Ar^+ spectra for the lowest perpendicular 
incident velocities which can be explained by a unique above-surface process involving the 
L-vacancy and two electrons from the resonantly populated shells. 

HCI beams have been deemed a candidate for future surface modification techniques for 
some time. It has been demonstrated that single ions can give rise to nanoscale size features 



on certain surfaces ||36|. Also sputter yields on insulators could be significantly enhanced by 
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using slow HCIs instead of fast singly charged projectiles. At very low kinetic energies, the 
energy deposition concentrates on a very small area which extends approximately one lattice 
constant in the vicinity of the first bulk layer. In this manner, an energy of several keV can 
be carried into this zone where it might be converted into activation energy for processes 
like sputtering, crystal growth and surface catalysis. Research in this field is under way and 
first results have been presented already. 
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FIGURES 

FIG. 1. Target chamber geometry, (a) The HCI beam colhdes with the target sm'face at an 
incident angle G. Electrons are emitted at an observation angle ^ with respect to the target 
surface towards the spectrometer entrance lens. In most experiments, we chose + ^ = 90°. For 
^ ^ 0°, electron emission originating from all interaction phases can be detected, (b) At grazing 
observation ^' = 0°, the detection of electrons emitted from below the surface is suppressed due to 
their damping on the long path through the bulk towards the spectrometer lens. Only electrons 
emitted from the above-surface zone are not affected such that spectral regions stemming from 
transitions of incoming or reflected projectiles gain intensity with respect to regions ascribed to 
below-surface emission. 

FIG. 2. Experimental spectra of Ar^+ projectiles impinging at a fixed incident angle of 
Q = 45° and various energies on a n-Si(lOO) surface. All spectra in this paper are normalized 
on the portrayed peak region between 160eV and 240eV. The error at the peak maximum of 
the (-Ej^in ~ ^'^^ eV)-spectrum is smaller than 0.3%. With increasing incident energy E^^^, the 
224eV-peak begins to outgrow the 232eV-peak and the main peak moves from 211eV to 208eV. 

FIG. 3. Ar^"*" spectra at E^^.^ = 121 eV and varying incident angles G. The intensities of the 
two subpeaks on the shoulder above 211eV alter with v± in the same way as in Fig. |. For the two 
geometries at ^' = 0°, a second broad peak region evolves around 198eV which we associate with 
emission from partially ionized incoming or reflected particles. 

FIG. 4. Ar^+ spectra at Ey^^^ = 9 eV and as parameter. The perpendicular velocity compo- 
nent v± is minimal in all four measurements and Gggp thus stays roughly constant due to the image 
potential. Due to our normalization method, the high-energy shoulder above 211eV is untouched 
by changing the nominal incident angle 0. At ^' = 0°, the broad region around 198eV pops up 
again. 
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FIG. 5. Series of Ar^+ spectra taken at fixed angles G = 90° and ^' = 0° with varying E^^^. 
For -Ejjjj^ = 9 eV and 121eV, we can still distinguish two separate peaks around 198eV and 211eV. 
For higher Ey^^^, the 211eV-peak disappears and the broad structure washes out while extending 
to lower energies. The 232eV-peak can only be found for minimum v±. 

FIG. 6. This histogram displays the energetic positions of LMM Auger transitions originat- 
ing from neutral "hollow" Ar^+ atoms comprising all possibles (2p^3s^p^d^)-configurations with 
hm = x + y + z<9. Angular coupling shifts have been neglected. The neutralization for con- 
figurations possessing less than q = 9 electrons in the M-shell has been established via {q — um) 
"spectator" electrons in the A^-shell. Each transition has been weighted by unity discarding tran- 
sition rates and statistical factors due to different subshell occupations. 

FIG. 7. Potential V^q (schematic) induced by an HCI core which is embedded into an electron 
gas. Vq^^'^ denotes a free ionic Coulomb potential. By adding outer spectator electrons to the atom, 
a screened Coulomb potential Vc'^^t"' forms which is implemented in our Cowan code calculations 
(spectator electron model). At large distances from the nucleus, the potential V^q seen by an active 
electron merges into the bottom of the valence band Vq. Electrons occupying states above Vq cannot 
be bound by the ionic core. They are lost to the valence band continuum. As level binding energies 
EJ}'' approach Vq, level shifts AE'^' increasingly deviate from the spectator electron model. 

FIG. 8. Binding energy ranges of hollow atom Ar^"^ 3£ subshells as a function of the total 
M-shell population. Within the spectator electron model, the 3d shell moves above the bottom of 
the valence band Vq as soon as more than five electrons populate this level. For comparison, we 
added the 2p level binding energies of C^^ and Ne^"^ as a function of the total L-population ni 
and a DFT calculation ^ for the Ne^+-Al system. 
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FIG. 9. Simulation of the experimental (cf. Fig. ^) Ar^+ and Ar^"*" spectra at = 45° for 
-^kin ~ ^ 121eV (b) and lOeV (c). For the {E-^[^ = 10 eV)-simulation, a statistically 

distributed M-shell population with (um) = 2 at the bulk entrance (see text) is implemented. The 
intensity ratios between the LMM transition subgroups as a function of E^^^ display the same 
systematics as the experiment. For the E-^^^^ = 2 keV, we enforce a correction for elastic and 
inelastic bulk damping by convolving the original spectrum with an exponential function of decay 
length of 3 a.u.. 

FIG. 10. Experimental spectra of metastable Ar^+ projectiles impinging at a fixed incident 
angle of G = 45° and varying E-^^^ on a n-Si(lOO) surface. Despite the initial 3s electron, the 
general shape of the spectra is very similar to the corresponding Ar^+ spectra in Fig. ^. Only the 
3ss- and 3sp region intensities are slightly enhanced. The small peak at 247eV occurs at minimal 



v± and can be assign to above-surface LMVw emission |35| 



FIG. 11. Experimental spectra of metastable Ar^+ projectiles impinging at E-^^^ = 8 eV and 
different incident angles on a n-Si(lOO) surface. All spectra exhibit the small peak at 247eV as 
expected for the minimum perpendicular projectile velocity components v±. At ^' = 0°, the broad 
peak around 198eV pops up again, cf. Figs. ^-|5|. In the (0 = 5°,^* = 0°)-geometry, its intensity 
is reduced with respect to the 211eV-peak as compared to the corresponding Ar^'^ spectrum in 
Fig.§ 
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TABLES 



process 


rfilled 


%n = 9 eV {{um) = 2) 


^kin = 121 eV 


%n = 2 keV 


3ss 


3.31 X 10^2 


0.8% (2.0|4.5|0.1) 


1.0% (2.0|4.3|0.1) 


2.0% (2.0|4.1|0.1) 


3sp 


5.29 X 10^3 


15.9% (1.6|5.1|0.0) 


17.1% (1.7|5.1|0.0) 


22.1% (2.0|5.0|0.0) 


3pp 


1.98 X 10^^ 


72.5% (1.4|5.5|0.0) 


70.0% (1.5|5.5|0.0) 


66.8% (2.0|5.4|0.0) 


3sd 


6.20 X 10^^ 


5.3% (1.2|0.5|1.4) 


7.2% (1.2|0.4|1.5) 


7.4% (2.0|0.3|1.4) 


3pd 


1.65 X 10^^ 


4.6% (0.6|1.5|1.4) 


3.7% (0.8|1.4|1.4) 


1.4% (1.5|1.2|1.2) 


3dd 


4.13 X lO^'' 


0.8% (0.4|0.4|2.4) 


1.0% (0.5|0.2|2.3) 


0.4% (1.2|0.1|2.2) 



TABLE L Monte Carlo simulation results on LMM processes for Ar^^ impinging on Si(lOO) 
with ^i^j^ = 9 eV, 121eV and 2 keV. T^j^f^ gives the LMM rate for a filled M-sheh as required for 
the implemented fit formula [^. For each simulation, we list the relative intensity and, in brackets, 
the average (n3s|n3p|n3(i)-configuration at the time of LMM decay which provides information 
about the evolution of the subsurface cascade. 



process 


^kin = 9 eV {{um) = 2) 


^kin = 121 eV 


^kin = 2 keV 


3sdd 


66.2% (0.2|0.4|2.4) 


81.8% (0.2|0.2|2.4) 


11.7% (1.0|0.1|2.2) 


3pdd 


16.3% (1.0|0.2|2.8) 


21.0% (1.0|0.2|1.5) 


17.9% (1.5|0.1|2.4) 



TABLE IL Monte Carlo simulation results on MMM processes for Ar^+ impinging on Si(lOO) 
with E-^^^ = 9 eV, 121eV and 2 keV. The table lists the average occurrence of each transition type 
and, in brackets, the average M-sublevel population at the time of MMM emission. Other sCK 
transitions are energetically forbidden. 
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pMCy r„-li 


Y.MCV r„-li 




-pMCV r„-li 
J- tot [s J 


n 
u 


Q Q9 V 1 nil 
y.yz X iu 


9 n7 V 1 ni2 

Z.Ui X iU 


(; CI s- ■] r\13 
D.Dl X iU 


Q 1 n V 1 ni4 

o.iU X iU 


1 

i 


1 91 V- 1 ni-^ 

i.Zi X J.U 


9 f- 1 ^ -1 ril3 
z.-J 1 X J.U 


1 Q V- 1 ni3 

V.Jo X lA) 


J . Uo X J.U 


9 




o.ZD X lU 


1 90 \/ 1 ni4 
i.ZZ X lU 


1 /1 1 N/ 1 ni5 

i.41 X iU 


Q 

o 




A Ad \,' 1 nl3 
4. 40 X iU 


1 /I /I 1 ni4 

i.44 X iU 


1 7n \/ 1 nl5 
i. ( U X iU 


4 




6.53 X 10^'' 




2.61 X 10^5 


5 




5.78 X 10^^ 




1.74 X 10^5 


6 




4.48 X 10^"^ 




9.17 X 10^"^ 


7 




3.27 X 10^^ 




3.27 X 10^^ 



TABLE III. MCV rates for the Ar^+/Si system. The table lists MCV transition rates per spin 
state rl^*^^ for each M-sublevel and the overall MCV rate T^f^ taking into account occupation 
statistics as evaluated by DFT calculations, um gives the initial number of M-electrons. The rates 
refer to initial M-shell ground state configurations. For um = 0, MCV processes filling the 3d 
level possess by far the highest rates. As the subsurface cascade proceeds and um > 4, the 3d level 
vanishes and the MCVs into the 3p level rapidly populate the M-shell. 
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